IPPP/06/25 
DCPT/06/50 
CERN-PH-TH /2006-070 
,hep-ph/0607104| 



Combined LHC/ILC analysis of a SUSY scenario with 

heavy sfermions 

K. Desch'^, J. Kalinowski^ G. Moortgat-Pick% K. Rolbiecki^ W.J. Stirling'^ 

°'Physikalisches Institut, Universitdt Bonn, D-53115 Bonn, Germany 
^Instytut Fizyki Teoretycznej, Uniwersytet Warszawski, PL- 00681 Warsaw, Poland 
'^TH Division, Physics Department, CERN, CH-1211 Geneva 23, Switzerland 
^IPPP, University of Durham, Durham DHl 3LE, UK 

February 2, 2008 

Abstract 

We discuss the potential of analyses at the Large Hadron Collider and the planned 
International Linear Collider to explore low-energy supersymmetry in a difficult region 
of the parameter space characterized by masses of the scalar SUSY particles around 
2 TeV. Precision analyses of cross sections for light chargino production and forward- 
backward asymmetries of decay leptons and hadrons at the first stage of the ILC with 
^/s = 500 GeV, together with mass information on X2 ^"^^ squarks from the LHC, allow 
us to determine the underlying fundamental gaugino/higgsino MSSM parameters and to 
constrain the masses of the heavy, kinematically inaccessible sparticles. No assumptions 
on a specific SUSY-breaking mechanism are imposed. For this analysis the complete spin 
correlations between production and decay processes must be taken into account. 



1 Introduction 

Supersymmetry (SUSY) is one of the best-motivated candidates for physics beyond the Stan- 
dard Model (SM). Low-energy SUSY is well- motivated since it stabilizes the electroweak scale, 
provides quantitatively accurate unification of gauge couplings as well as a promising cold-dark- 
matter candidate. Moreover, electroweak precision data [1] and cosmology bounds [2] seem 
to indicate [3] that at least some of the electroweakly interacting SUSY particles should be 
rather light and accessible at future colliders. However, since the mechanism of SUSY breaking 
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is unknown, supersymmetric extensions of the Standard Model contain a large number of un- 
known parameters, e.g. 105 in the Minimal Supersymmetric Standard Model (MSSM). Specific 
assumptions on the SUSY-breaking mechanism, in particular about the unification of parame- 
ters at the grand-unification (GUT) scale, considerably reduce the number of free parameters, 
e.g. in the constrained MSSM, often referred to as mSUGRA, where we end up with only 
four new parameters (and one sign) specified at the unification scale. Experiments at future 
accelerators, the Large Hadron Collider (LHC) and the International Linear Collider (ILC), 
will have, however, not only to discover SUSY but also to determine precisely the underlying 
SUSY-breaking scenario with as few theoretical prejudices as possible. 

Particularly challenging are scenarios where the scalar SUSY particle sector is very heavy 
as required, for instance, in focus-point scenarios (FP) [4] where the gaugino masses are kept 
relatively small while squarks and sleptons might be too heavy for a direct observation at 
the ILC. It is therefore particularly interesting to verify whether the interplay of a combined 
LHC/ILC analysis [5] could shed hght on models with heavy sfermions. 

Many methods have been worked out to derive the SUSY parameters at collider experi- 
ments [6-10]. In [11-16] the chargino and neutralino sectors have been exploited at the ILC to 
determine the MSSM parameters. However, in most cases only the production processes have 
been studied and, furthermore, it has been assumed that the masses of the virtual scalar par- 
ticles are already known. In the case of heavy scalars such assumptions cannot be applied and 
further observables have to be used to determine the underlying parameters. Studies have been 
made to exploit the whole production-and-decay process, and angular and energy distributions 
of the decay products in chargino as well as neutralino channels have been studied in [17-20]. 
Since such observables depend strongly on the polarization of the decaying particle, the spin 
correlations between production and decay can have a large influence and have to be taken 
into account. Exploiting such spin effects, it has been shown in [21] that, once the chargino 
parameters are known, useful indirect bounds for the mass of the heavy virtual particles could 
be derived from forward-backward asymmetries of the final Icpton Afb(^)- 

In this paper we discuss a FP-inspired scenario characterized by a ~ 2 TeV scalar particles 
sector. In addition, the neutralino sector turns out to have very low production cross sections 
in e^e" collisions, so that it might not be fully exploitable. Only the chargino pair production 
process has high rates at the ILC and all information obtainable from this sector has to be 
used. In order to assess the possibility of unraveling such a challenging new physics scenario, 
our analysis is performed entirely at the EW scale, without any reference to the underlying 
SUSY-breaking mechanism. We measure at the LHC and at the ILC with ^/s = 500 GeV 
the masses, cross sections and spin-dependent forward-backward asymmetries and analyze the 
potential of a multiparameter fit to determine the underlying parameters. 

The paper is organized as follows. In the next section we first present the studied pro- 
cess, chargino production with leptonic and hadronic decays. We briefiy introduce the spin 
formalism, which is needed for the evaluation of spin-dependent observables. In section 3 the 
FP-inspired scenario is defined and the expected experimental results at the LHC and the ILC 
are discussed. In section 4 we perform our numerical analysis and determine the SUSY param- 
eters based on the experimental input. An attempt at testing the SU{2) symmetry relation for 
the selectron and sncutrino masses using the available information on the squark masses from 
the LHC and the forward-backward asymmetry measured at the ILC in hadronic decay modes 
is also discussed. Section 5 summarizes the results. 
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e+ xt e+ xi e+ xi 




Figure 1: Feynman diagrams for production and for leptonic and hadronic decays of charginos. 

2 Strategy overview 

2.1 Chargino and neutralino sector 

We study chargino production 

e- + e^^xt + Xi, (1) 
with subsequent leptonic and hadronic decays 

xt Xi + ^^ + i^ and Xi + qa + qu, (2) 

Xi Xi + ^~' + and Xi + C[d + Qu, (3) 

where i = e, fi, Qu = u, c, = d, s. The corresponding Feynman diagrams are shown in Fig. [TJ 
The production process contains contributions from 7- and Z'^-exchange in the s-channel and 
from //-exchange in the t-channel. The decay processes are mediated by W^, ii,, z> or by qdL, QuL 
exchange; contributions from Higgs boson exchanges to the production and decay are neghgibly 
small for the first and second generation fermions. For notations, couplings and conventions 
see, for instance [22]. 

The neutralino mass eigenstates are defined as = ^ij^ji where Nij are the elements of 
the unitary 4x4 matrix which diagonalizes the neutral gaugino-higgsino mass matrix in the 
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basis X^{B^,W^,H^, Hi): 
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mz sin /? sin \ 
—mz sin /3 cos 9w 



-1^ 




(4) 



where mz denotes the mass of the Z° boson, Mi, M2 are the SU{2) gaugino mass 

parameters, /i is the Higgs mass parameter and tan/3 = V2/V1, where Vi^2 are the vacuum 

expectation values of the two neutral Higgs fields. The chargino mass eigenstates xt — (--) 

are defined by xt = VnW^ + ^2^"*" and xj = Ujiw~ + Uj2h~ . Here w"^ and h"^ are the two- 
component spinor fields of the wino and the charged higgsinos, respectively. Furthermore, Uij 
and Vij are the elements of the unitary 2x2 matrices, which diagonalize the chargino mass 
matrix: 

M2 mwV^ sin P 
mwV^ cos P jj, 

where mw denotes the mass of the bosons. 



M 



c 



(5) 



2.2 Spin formalism 

We study the production process including the subsequent leptonic and hadronic decays. Since 
in our scenario charginos are very narrow (r^± = 2.3 keV), the narrow- width approximation is 
appropriate and contributions from off-shell channels are negligible; this approximation can be 
tested with several Monte Carlo event generators that include off-shell effects as well as spin 
correlations [23] . The process can therefore be split into the chargino production and the decay 
processes. In order to exploit spin-dependent observables, e.g. forward-backward asymmetries 
of the final leptons and quarks, however, the full spin information of the decaying charginos has 
to be taken into account. In the following we briefiy summarize the required spin formalism. 
More details as well as the explicit analytic expressions for the chargino production spin density 
matrix and the decay processes can be found in [17, 18]. 
The amplitude for the whole process is 

T = A(x+)A(x7) T"p''Tn,,Tn,x,, (6) 

Xi,\j 

with the hehcity amphtude for the production process Tp^^ and those for the decay processes 
TD,Xi, TD,Xj, and the propagators A(x^) = l/[si ~ mf + imiVi]. Here Aj, Sj, m^, Ti denote the 
helicity, four- momentum squared, mass and width of xf. The amplitude squared 

|rp = |A(x+)r|A(x-)rp''''''^'^'^Py.A.Py,A, (sum convention used) (7) 

is thus composed of the (unnormalized) spin density production matrix 

pP,XiX,xJ^, = T^^'^Tp^'^* (8) 



4 



of xtj^ and the decay matrices 

Px'^K = TD,xJh,x'^ and p^,^^^ = TD,x,T^,xy (9) 

Interference terms between various hehcity amphtudes preclude factorization in a production 
factor XIajAj I^p' times a decay factor J2xiXj\'^D,XiXj\'^- 

The spin density production matrix p^'^^^^^i^j can be decomposed into four parts (for details 
see [17]): 

where P denotes a contribution which is independent of chargino polarizations, (Sf,) depends 
on the polarization of one of the charginos, and T^ab on both; a,b = 1,2,3 denote the components 
(transverse and longitudinal) of the spin vectors. Likewise, the decay matrices p?,._ and p?, 
can each be separated into two parts as 



PX'X, = Dy^X.+Kyx.- (11) 

The amplitude squared |Tp of the combined processes of production and decays can schemat- 
ically be written as (with helicity indices suppressed): 

|Tp ~ PD,D, + Sf + Sf A + Sr,Sf , (12) 

The first product in Eq. f|T2l) is independent of spin correlations between production and decay. 
The second and third terms describe the correlations between the production and the decay 
process either of or xj decay and, in the last term correlations between both decay processes 
are included. 

• In the first term of Eq. f|T2|) only scalar products appear, which can be expressed by the 
Mandelstam variables s, t, u for the production and decay processes. This is the only term 
that survives in the total cross section, i.e. when all the angles are integrated over. 

• In the second (third) term of Eq. f|T2l) the spin vectors relate quantities from the pro- 
duction with those from the decay process. These scalar products cannot be expressed 
by Mandelstam variables. They contain the angle between the incoming electron and 
the outgoing lepton or quark in the laboratory system. These terms contribute to spin- 
dependent observables as, for example, to the forward-backward asymmetry of the final 
leptons and quarks. 

• The last term of Eq. (fT2|) involves the spin vectors of both charginos and leads to spin 
correlations between the two decay chains of xt and x] ■ 

If the decay of only one chargino, e.g. xt , is considered, one has to sum over the spin of x] so 
that in Eq. ^ Dj = l and = 0. 
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2.3 Strategy for the determination of the SUSY parameters 



Our aim is to demonstrate the power of forward-backward asymmetries in determining the 
masses of kinematically inaccessible heavy sleptons. Note, however, that the leptonic forward- 
backward asymmetry involves the masses of both sneutrinos and selectrons. Therefore the 
analysis is performed in steps. First we exploit the chargino and neutralino masses and the 
chargino production cross sections times decay branching fractions, to constrain fundamental 
parameters of the chargino and neutralino sectors and the sneutrino mass. We then show how 
the obtained limits on the sneutrino mass and other parameters can be improved by employing 
the forward-backward asymmetries, measured in leptonic chargino decays at the ILC. This, 
however, requires the assumption on the SU{2) mass relation for slepton masses. In the last 
step we make an attempt at testing the SU{2) symmetry relation for the selectron and sneutrino, 
using the available information on the squark masses from the LHC, and include in addition 
the forward-backward asymmetry measured at the ILC in hadronic decay modes. 

3 Case study with heavy sfermions 
3.1 Parameters of the chosen scenario 

The following mSUGRA parameters, taken at the GUT scale except for tan j3, define the MSSM 
scenario: 

mi/2 = 144 GeV, mo = 2 TeV, = GeV, tan/? = 20, sgn(^) = +1. (13) 

However, our analysis is performed entirely within the general MSSM framework, without any 
reference to the underlying SUSY-breaking mechanism. The parameters at the EW scale are 
obtained with the help of the SPheno code [25] for rrit = 178 GeV; furthermore it has been 
checked with the code micrOMEGA [26] that the lightest neutralino provides a relic cold- 
dark-matter density consistent with cosmological data. The low-scale gaugino/higgsino/gluino 
masses as well as the derived masses of SUSY particles are listed in Tables [1] and [2l The 
charginos and neutralinos as well as the gluino are rather light, whereas the scalar SUSY 
particles have masses about 2 TeV. 



Ml 


Ms 


Ms 


/i 


tan/? 








X2 








60 


121 


322 


540 


20 


117 


552 


59 


117 


545 


550 


416 



Table 1: Low-scale gaugino/higgsino/ta.n/3 MSSM parameters, and the resulting chargino, neu- 
tralino and gluino masses (all masses are given in GeV). 



3.2 Expectations at the LHC 

As can be seen from Tables [T] and [21 all squarks are kinematically accessible at the LHC. The 
largest squark production cross section is for ti 2. However, with stops decaying mainly to gt 
(with BR{ii^2 gt) ~ 66%), where background from top production will be large, no new 



6 



rrih 


mH,A 


mH± 




















119 


1934 


1935 


1994 


1996 


1998 


1930 


1963 


2002 


2008 


1093 


1584 



Table 2: Masses of the SUSY Higgs particles and scalar SUSY particles (all masses are given 
in GeV). 



interesting channels are open in their decays. The other squarks decay mainly via gq, but 
since the squarks are very heavy, mg^R ~ 2 TeV, precise mass reconstruction will be difficult. 
Therefore we conservatively assume that the squark masses can be measured with an error of 50 
GeV. Our results do not depend sensitively on this assumption since the mere indication that 
the scalar quarks are very heavy will be sufficient for narrowing the experimental uncertainty 
on the slepton sector from the ILC measurements. 

In this scenario, the inclusive discovery of SUSY at the LHC is possible mainly because 
of the large gluino production cross section. Therefore several gluino decay channels can be 
exploited. The largest branching ratio for the gluino decay in our scenario is a three-body decay 
into neutralinos, BR{g X2^&) ~ 14%, with a subsequent three-body leptonic neutralino decay 
BR{x2 ~* ^ = e, /i of about 6%, see Table [3l In this channel the dilepton edge will be 

clearly visible, since this process has low backgrounds [5] . The mass difference between the two 
light neutralinos can be measured from the dilepton edge with an uncertainty of about [24] : 

5(m^o - m^o) ~ 0.5 GeV. (14) 



Mode 


g xlhh 


9 ^Xi QuQci 


Xt XiQaQu 


xl ^ xT^- 


ti,2 gt 


Xi ^Xi^ 


BR 


14.4% 


10.8% 


33.5% 


3.0% 


66% 


11.0% 



Table 3: Branching ratios for some important decay modes in the studied MSSM scenario, 
i = e, fi, T, qu = u,c, qd = d,s. Numbers are given for each family separately. 



Other frequent gluino decays are into the light chargino and jets, with BR{g Xi ll') ~ 20% 
for qq' in the ffist two families, and about 3% in the third, with a subsequent leptonic chargino 
decay BR{xi XA^^^t)i^ = e, of about 11%. However, exploiting this channel for the 
chargino-neutralino mass difference measurement is very difficult. First, because of the escaping 
neutrino, and second, because of a genuine 3-body chargino decay in our scenario. To our 
knowledge, the only attempt to determine chargino mass at the LHC required —>■ Xi'^^ ~^ 
Xi^'^^e with on-shell arriving at a statistical accuracy of ~25 GeV [27]. 

In both gluino decay channels the spin measurements via angular correlations [19, 20] of 
decay products should provide evidence for the spin 1/2 character of the intermediate particles 
assuring us that the underlying SUSY scenario is realized. 

Finally, the gluino mass can be reconstructed in a manner similar to the one proposed 
in [29], where the SPSla scenario is analyzed. Although our scenario is different, the precison 
in both is limited by systematic uncertainties due to hadronic energy scale and a similar relative 
uncertainty of ~2 % can be expected. 
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3.3 Expectations at the ILC 

At the ILC with ^/s < 500 GeV, only hght charginos and neutrahnos are kinematically ac- 
cessible. However, in this scenario the light neutralino sector is characterized by very low 
production cross sections. For example, at 500 GeV and {Pe-,Pe+) = (—90%, +60%) beam 
polarization we obtain (j{xiX2) = 0-93 fb, cr(x2X2) = 0-49 fb; for other beam polarization 
and/or lower collider energy the cross sections are even lower. This is due to the almost pure 
gaugino nature of light neutralinos (x? ~ 99%i?°, X2 ~ 97%Vr°) with suppressed couplings to 
the s-channel Z-boson, while the t- and w-channel selectron exchange is small because of the 
heavy selectron mass. Only the X3X4 channel (because of opposite CP parities of Xs and xt) 
could have an appreciable cross section above its threshold i/s ~ 1100 GeV. 

Only the chargino pair production process has high rates at the ILC and all information 
obtainable from this sector has to be used. We constrain our analysis to the first stage of the 
ILC with a/s < 500 GeV and study only the xfxi production 

e+e" xtxl, (15) 

with subsequent chargino decays 

Xi X?e"z/e, xldu and Xi^c (16) 

and the corresponding charge conjugate xf decays, for which the analytical formulae, including 
the complete spin correlations, are given in a compact form, e.g. in [17]. 

3.3.1 Mass measurements 

The chargino mass can be measured at y/s = 350 and 500 GeV in the continuum, with an error 
of about 0.5 GeV [28,30]. This can serve to optimize the ILC scan at the threshold [31] which, 
because of the steepness of the s-wave excitation curve in xfxi production, can be used to 
determine the light chargino mass very precisely, to about [30]: 

m-± = 117.1 ±0.1 GeV. (17) 

The mass of the lightest neutralino m^o can be derived via the decays of the light chargino, 
either from the energy distribution of the lepton i~ {BR{xi ^ Xi^^^i) ~ 11%, see Table [3]) or 
from the invariant mass distribution of the two jets in hadronic decays {BR{x^ —>■ XiQdQu) ~ 
33%, see Table [3]). We take [30] 

m^o = 59.2 ± 0.2 GeV. (18) 

Together with the information from the LHC, Eq. (fT4|) . a mass uncertainty for the second 
lightest neutralino of about 

m^o = 117.1 ±0.5 GeV (19) 

can be assumed. 
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3.3.2 Experimental uncertainties for the cross sections 

Table H] lists the expected chargino production cross sections and forward-backward asym- 
metries for different beam energies and polarization configurations, derived with the nomi- 
nal values of parameters. Experimentally we identify the chargino pair production process, 
e+e" — *• XiXi: ill the fully leptonic {i~^i^Xi^~^Xi) semileptonic {ii^Xiqq'Xi) filial states 
(where £ = e,/i). We estimate an overall selection efficiency of 50%. For both final states, 
W^W~ production is expected to be the dominant SM background. The fully leptonic chan- 
nel is more challenging, owing to the absence of mass constraints. Its efficient selection needs 
further experimental study. However, the fully leptonic channel is not essential for this anal- 
ysis, as its relative contribution to the overall rate is only ~ 14%. For the semileptonic 
(sic) final state, this background can be efficiently reduced from the reconstruction of the 
hadronic invariant mass. In Table HI we list cross sections multiplied by the branching fraction 
Bsic = 2 X BR{xl ^ Xlmu) X BR{xi ^ xT^) + [BR{Xi ^ xT^)f ~ 0.34 (£, g,, g„ 
include the first two families of leptons and quarks) including an Cgic = 50% selection efficiency. 
The error includes, added in quadrature, the statistical uncertainty based on numbers of iden- 
tified events for C = 200 fb~^ in each polarization configuration, {Pe-,Pe+) = (— 90%, -|-60%) 
and (+90%, —60%), and a relative uncertainty in the polarization of APe±/Pe± = 0.5% [33]. 
Since the production rates are high, the total uncertainties are rather small; see Table |H 



^/GeV iPe-,Pe+) 


^{xtxl)/Po 


(^{xtxl)Bsic e.ic/fb 


AMn/% 


Afb(c)/% 


350 (-90%, +60%) 


6195.5 


1062.5+4.0 


4.42+0.29 


4.18+0.74 


(+90%, -60%) 


85.0 


14.6+0.7 






500 (-90%, +60%) 


3041.5 


521.6+2.3 


4.62+0.41 


4.48+1.05 


(+90%, -60%) 


40.3 


6.9+0.4 







Table 4: Cross sections for the process e^e^ — > XiXi ^^'^ forward-backward asymmetries 
(A-fb) in the leptonic Xi ~^ Xi^~^ ^nd hadronic Xi ~^ x\sc decay modes, for different beam 
polarization P^- , Pe+ configurations at centre-of-mass energies i/i = 350 GeV and 500 GeV at 
the ILC. Concerning the errors, see text. 



3.3.3 Experimental uncertainties for the forward— backward asymmetries 

Figure [2] shows the expected forward-backward asymmetries measured in the leptonic and 
hadronic decay channels as functions of the sneutrino mass; for illustration, the dashed line 
in the left panel demonstrates that spin correlations between production and decay must be 
taken into account for a proper interpretation of the experimental data. In the case of leptonic 
decay the SU (2) relation on slepton masses 

"^gL = + ^1 cos(2/5)(-l + sin^ 9^) (20) 

has been assumed, while for the hadronic decay the squark masses are taken to be fixed by the 
LHC measurement. 
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The forward-backward asymmetry is defined as: 

APB = (21) 

where ap and ctb are the acceptance-corrected cross sections. The statistical error on A-p-g, 
assuming a binomial distribution, is given by 

A(Afb)^*"* = 2v/e(l - e)/N, (22) 

where e = (7f/(o"f + o"b) and N denotes the number of selected events. Errors due to uncer- 
tainties of beam polarization and squark masses are negligible. 

For the leptonic forward-backward asymmetry, a-p and ctb are the acceptance-corrected 
cross sections for the e~/fi~ from Xi being observed in the forward (F) and backward (B) 
hemisphere, respectively, for i'i^XiQQ'Xi and i~uXi^~^^Xi final states. For the ^"^i^XiQQ'Xi final 
state, in which comes from xt decay, the condition is that the positive lepton is observed 
in the backward (F) or forward (B) direction. 

The hadronic forward-backward asymmetry has been analyzed if one fiavour of the hadronic 
2-jet system can be identified. Since the expected vertex-detector performance is excellent 
at the ILC, charm-tagging from secondary vertices and displaced tracks can be envisaged. 
This allows us to measure the forward-backward asymmetry of the c/c-jet in each of the 
different semileptonic decays i'^i'XiCSXi and i'uXiCSXi- Note that the charge sign of the c-jet 
can be inferred from the lepton charge of the other chargino decay and does not have to be 
measured through jet- or vertex-charge techniques. Taking the results from [34], c-jets can 
be identified with an efficiency of 50% at a purity of 80% in Z-decays at rest. The purity 
from chargino and W decays is larger, since the ratio of true charm/non-charm jets is ~ 1/3 
(compared to ~ 1/5 for Z decays). In this analysis, we assume an overall selection efficiency 
for iuXiCSXi of 20%, corresponding to a c-tag efficiency of 40% for a selection efficiency of 50%. 
In Table m the asymmetries are listed only for the (Pe-,Pe+) = (—90%, +60%) case, since the 
cross sections for the opposite polarization are very small and the statistical errors become very 
large. Consequently we do not include them in the following analysis. 



4 Parameter determination 

In the following we apply multiparameter fits to determine the underlying SUSY parameters. 

• In the first step we use only the masses xf, Xi? X2 the chargino pair production 
cross section, including the full leptonic and the semileptonic decays as observables. A 
four-parameter fit for the parameters Mi, M2, fi and m,> has been applied. 

• In the second step we include as an additional observable the leptonic forward-backward 
asymmetry. Only the semileptonic and purely leptonic decays were used. The SU{2) 
relation between the two virtual masses mc, and m^^ has been applied as an external 
constraint. 

• As an attempt to test the SU{2) mass relation for the slepton and sneutrino masses, in 
the last step both the leptonic and hadronic forward-backward asymmetries have been 
used. A six-parameter fit for the parameters Mi, M2, /i, mj>, and tan (3 has been 
applied. 
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Figure 2: Forward-backward asymmetry of e~ in the process e+e" — ^ XiXi > Xi ~^ Xi^~^ (^^ft 
panel) and of c- jet in the process e+e" xtXi > Xi ~^ Xi-^c (right panel), shown as a function 
ofmc,^ at ^/s = 350 GeV and polarized beams {Pe-, Pe+) = (—90%, +60%). In the left panel the 
mass of the other scalar virtual particle, rrie^, which contributes to the decay process, has been 
assumed to fulfill the SU (2) mass relation, Eq. ^20). while in the right panel the mass of the 
squark is kept fixed as measured at the LHC. For nominal value ofrriy^ = 1994 GeV the expected 
experimental errors are shown, see Eq. For illustration only, the dashed line in the left 
panel shows that neglecting spin correlations would lead to a completely wrong interpretation of 
the experimental data. 



4.1 Parameter fit without using the forward— backward asymmetry 



We use as observables the masses m-+ 



m^o, m^o and the polarized chargino cross section 
multiphed by the branching ratios of semileptonic chargino decays: a{e^e~ —>■ XiXi) ^ Bsic, 
with Bsic = 2 X BRixt ^ X^qaqu) x BRixi ^ X?^-^>) + [BRixi x'l^-i^)? ~ 0.34, £ = e,/i, 
Qu = u, c, Qd = d, s, with selection efficiency 50%, as given in Table HI Note that the chargino 
branching ratios are not sensitive functions of sfermion masses mediating their decays, since we 
know from the LHC that sfermions are very heavy. We take into account the la statistical error 
based on £ = 200 fb^^ for each polarization configuration, a relative uncertainty in polarization 
of APf>±/Pe± = 0.5% [33] and an experimental efficiency of 50%; see Table HI 

We apply a four-parameter ffi for the parameters Mi, M2, /i and rric,^ for fixed values of 
tan/3 = 5, 10, 15, 20, 25, 30, 50 and 100. Fixing tan (3 is necessary for a proper convergence of 
the fitting procedure [32] because of strong correlations among parameters. We perform a 
test 

a{ij) - cr{ijy^^ 



X 



E 



j=LR.RL 
j=350,500 



rrii 



m 



thx 2 



(23) 



where in the ffist term i = LR, RL denote {e~,e~^) polarization configurations (—90%, +60%) 
and (+90%, —60%) respectively, and j = 350,500 denote the cm. energy in GeV; <j{ij), mi 
and AcT(zj), Amj are the corresponding cross sections and masses and their uncertainties; see 
Tableland Eqs. (fT71)-(fT9|). 

It turns out that for tan/3 < 1.7 the measurements are inconsistent with theoretical predic- 
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Figure 3: Migration of la contours for tan (3 = 5, 10, 20, 50 with the other two parameters fixed 
at the values determined by the minimum of for each tan/3. 



tions at least at the la level. The corresponding la constraints from cross sections and light 
neutralino and chargino mass measurements for the underlying parameters are as follows 

59.4 < Ml < 62.2 GeV, 118.7 < M2 < 127.5 GeV, 

450 < ;u < 750 GeV, 1800 < m^^ < 2210 GeV. (24) 

Owing to the strong gaugino component of and X125 parameters Mi and M2 are deter- 
mined reasonably well, with a relative uncertainty of ~ 3% and ~ 5%. The higgsino parameter 
jj, as well as m^^ are determined to a lesser degree of precision, with relative errors of ~ 30% and 
10%. Note, however, that large errors are partly due to migration of the fitted central values 
of the parameters with tan/3. Figure [3] shows the migration of la contour^ in my^-M2 (left), 
M2-/i (middle) and M1-M2 (right), the other two parameters being fixed at the values deter- 
mined by the minimum of fo^' tan/3 changing from 5 to 10, 20 and 50. Beyond tan/3 = 50, 
the migration is negligible. Varying tan/? between 5 and 50 leads to a shift ~ 1 GeV of the 
fitted central Mi value and ~ 3.5 GeV of M2, effectively increasing their experimental errors, 
while the migration effect for and nii,^ is much weaker. 



4.2 Parameter fit including the leptonic forward— backward 
asymmetry 

We now extend the fit by using as additional observable the leptonic forward-backward asymme- 
try for polarized beams (—90%, +60%). We include final-state electrons and muons, assuming 
equal masses of selectrons and smuons, and we include decays of both charginos. The SU{2) 
relation between selectron and sneutrino masses has been assumed, see Eq. f l20l) . The param- 
eter ranges found in the previous step are scanned and accepted if xipg ^ 1 after inclusion of 
forward-backward asymmetry according to 

i 

^Note that these plots are 2-dim cuts of a 4-dim hypersurface for each tan /3 value and they may give a false 
impression that errors are smaller than that in Eq. (j24|) . 
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Figure 4: Two 2-dim cuts in {M1-M2} and {^-niy} planes of the la hypersurface before (dashes) 
and after (solid) including the leptonic ApB in the fit for fixed ta.n /3 = 20 and the other param- 
eters taken at the values determined by the minimum ofx^- 

where is defined as in Eq. fl2^ , and the sum over i includes ApB measured for both electrons 
and muons at cm. energies of 350 and 500 GeV. The terms Afb(^) and AApB^i) are the cor- 
responding experimental forward-backward asymmetries and their uncertainties; see Table |H 
For the forward-backward asymmetries the errors due to the uncertainty of beam polarization, 
although very small with respect to the statistical one, are also included in the test. The 
effect of including the leptonic FB asymmetry in the fit is illustrated in Fig. H] for two 2-dim 
cuts of the la hypersurface at fixed tan /? = 20. As expected, the range of rri;^ is most affected, 
although changes of contours are not large. However, the main virtue of including Aps, not 
visible in this Figure, is constraining tan (3. No assumption on tan (3 has to be made in the fit 
since for too small or too large a value of tan (3 the wrong value of Ap^ is predicted. As a result, 
including the forward-backward asymmetries in the multiparameter fit strongly improves the 
results. We find 

59.7 < Afi < 60.35 GeV, 119.9 < M2 < 122.0 GeV, 500 < < 610 GeV, 

1900 < m^^ < 2100 GeV, 14 < tan/? < 31. (26) 

In particular, tan (3 is constrained from below rather well. The constraints for the mass m^^ are 
improved by a factor of about 2 and for gaugino mass parameters Mi and M2 by a factor of 
about 5, as compared to the results of Section H?T] with unconstrained tan/3. The error for the 
higgsino mass parameter /i also decreases significantly. From Eq. (I26l) we obtain the following 
predictions for the heavy charginos/neutralinos 

506 < m^o < 615 GeV, 512 < m^o < 619 GeV, 514 < m-± < 621 GeV. (27) 

These sparticles are only kinematically accessible with non-negligible cross sections at a phase-2 
ILC. 
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4.3 Parameter fit including the hadronic and leptonic 
forward— backward asymmetries: test of SU{2) 

4.3.1 Parameter fit including the leptonic Ape 

In principle the SU (2) relation can be tested by employing the forward-backward asymmetries 
measured in hadronic and leptonic decay modes of produced chargino. With the constraints 
for the squark masses from the LHC, the hadronic forward-backward asymmetry could be used 
to control the sneutrino mass and the leptonic forward-backward asymmetry to derive con- 
straints on the selectron mass. However, with the foreseen experimental accuracies, testing the 
SU (2) relation turns out to be very challenging because our measurements do not sufficiently 
constrain all 6 parameters simultaneously. Therefore, we perform a scan of the parameter 
space and calculate Xa according to Eq. fl25|) . i.e. taking light chargino and neutralino masses 
(Eqs. fll7til9p ). and 4 cross section measurements and 2 leptonic forward-backward asymmetry 
measurements (entries in columns 3 and 4 of Table H]). From xipB = 1 we derive the following 
constraints: 

59.30 < Ml < 60.80 GeV, 117.8 < Ms < 124.2 GeV, 420 < < 950 GeV, 
1860 < nii,^ < 2200 GeV, mg^ > 1400 GeV, tan/5 > 11. (28) 

As we can see, without the SU{2) relation the upper limits on the selectron mass and tan/3 
cannot be established since a change of these parameters for high values can be compensated 
by small changes of other parameters. Limits for the parameter /z are also very poor. The 
parameters Mi and M2 are nevertheless quite well determined, with an accuracy of the order 
of a few per cent, thanks to tight experimental mass constraints on the light chargino and 
neutralinos. 

4.3.2 Parameter fit including the leptonic and hadronic ApB 

Including hadronic forward-backward asymmetry (two entries in the last column of Table H]) 
improves the constraints as follows: 

59.45 < Ml < 60.80 GeV, 118.6 < M2 < 124.2 GeV, 420 < < 770 GeV, 
1900 < < 2120 GeV, mg^^ > 1500 GeV, 11 < tan/? < 60. (29) 

The most significant change is for the sneutrino mass, for which error bars become smaller by ~ 
50%. Also an upper limit on tan /? is found, which has the effect of improving the upper limit on 
// significantly and slightly lowering the limits for Mi and M2. However we do not get an upper 
limit for the selectron mass. Nevertheless, the results for the selectron and sneutrino masses 
are consistent with the SU{2) relation. The hadronic forward-backward asymmetry would be 
much more useful with more precise measurements, which is very challenging experimentally. 

5 Conclusions 

We showed a method for determining the MSSM parameters in scenarios with heavy scalar 
particles where only a small part of the particle spectrum is kinematically accessible at the 
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ILC. Such scenarios appear very challenging, since only a very limited amount of experimental 
information is accessible about the SUSY sector. However, a careful exploitation of data leads 
to significant constraints for unknown parameters. A very powerful tool in this kind of analysis 
turns out to be the forward-backward asymmetry. The proper treatment of spin correlations 
between the production and the decay is necessary in this context. This asymmetry is strongly 
dependent on the mass of the exchanged heavy particle. If the SU{2) constraint is applied, the 
slepton masses can be determined to a precision of about 5% for masses around 2 TeV at the 
ILC running at 500 GeV, i.e. one eighth of the energy necessary for slepton pair production. 
Also the derived constraints on heavy chargino/neutralinos may provide the physics argument 
for a second stage of the ILC. 

The SU{2) assumption on the left-selectron and sneutrino masses could in principle be 
tested by combining the leptonic forward-backward asymmetry with the hadronic forward- 
backward asymmetry for which the squark masses are measured at the LHC. However, with 
current estimates for the efficiency and purity of charm tagging at the ILC, this test is not very 
stringent. With significantly better charm-tagging performance more sensitive tests could be 
performed. 

Our analysis stresses the important role of the LHC/ILC interplay, since neither of these 
colliders alone can provide us with the data needed to determine the SUSY parameters in such 
scenarios with heavy sfermions without tight model assumptions. 
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